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I. INTRODUCTION 



The measurement of longitudinal spin asymmetries in production of W ± bosons at the 
Relativistic Heavy Ion Collider (RHIC) will provide an essential probe for spin-dependent 
quark distributions at high momentum scales Q 2 [1]. At app center-of-mass energy y/s = 500 
GeV, about 1.3 x 10 6 VF + and W~ bosons will be produced by the time the integrated 
luminosity reaches 800 pb _1 . Due to the parity violation in the qqW coupling, this process 
permits non- vanishing single-spin asymmetries Al(£), defined for any kinematical variable 
f as 

da(p^p->WX) _ dajp^p^WX) 

A (c\ — d € (i\ 

d£ i d£ 

The Born-level expression for the asymmetry A L (y w ) with respect to the rapidity yw of 
the W boson is particularly simple if the absolute value of yw is large. In that case, the 
Born-level A^yw) reduces to the ratio Aq(x)/q(x) of the polarized and unpolarized parton 
distribution functions [2-4]. Furthermore, A^yw) tests the flavor dependence of quark 
polarizations. 

The original method for extracting the spin-dependent quark distributions out of the 
RHIC W boson data is based on the direct reconstruction of the asymmetry A L (y w ) [1, 5]. 
Unfortunately, such reconstruction is obstructed by specifics of the detection of W ± bosons 
at RHIC. First, neither Phenix [6] nor Star [7] detector at RHIC is hermetic. Therefore, 
it is not possible at those detectors to monitor the energy balance in particle reactions, 
even in the transverse direction with respect to the beams. Due to the lack of information 
about the missing energy carried by the neutrino, the determination of yw is in general 
ambiguous and depends on the assumptions about the dynamics of the process. Second, 
due to the correlation between the spins of the initial-state quarks and final-state leptons, 
the measured value of Ai(yw) is strongly sensitive to the experimental cuts imposed on the 
observed charged lepton. Hence, unless a theory calculation exists to reliably predict the 
distributions of the leptons from the W boson decay, the spin-dependent quark distribution 
functions cannot be determined with acceptable accuracy. 

In weak boson production at the Fermilab Tevatron pp collider, the mass Mw of 
the W boson has been measured to a great precision [8, 9]. The information on M w 
is extracted out of the lepton-level distributions da/rriT, dcr/dpTe, or da/dpr u , where 
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PTe = | Pre I an( ^ Vtv = |P7v| are the transverse momenta of the electron and neutrino, 
and rriT = \/2px e PTu — 2pr e • Ptv is the transverse mass of the lepton pair. The value of 
Mw can be determined from the analysis of the position and shape of the Jacobian peaks 
in each of these three distributions. For most of the events contributing near the Jacobian 
peak, the transverse momentum qr of the W boson is very small, of order of a few GeV. 
Therefore, at each order of the perturbative QCD calculation, there are large logarithmic 
contributions of the form o. s q^ 2 ln m (M^/g|,), where n = 1, 2, . . . and m — 0, 1, . . . , 2n — 1. 
These large logarithmic contributions can be summed to all orders in as by applying the 
transverse momentum resummation formalism formulated in the impact parameter space by 
Collins, Soper, and Sterman (CSS) [10]. Without such a resummation calculation performed 
at the decay lepton level [11], it would not be possible to determine Mw from the experi- 
mental data. This is because at hadron colliders the longitudinal momentum of the neutrino 
from the W boson decay is not observed by the detectors; hence, the invariant mass of the W 
boson cannot be measured directly. Instead, the determination of the W boson mass relies 
on the analysis of the shape of the Jacobian peaks in the lepton-level distributions sensitive 
to qr- Again, without a resummation calculation at the secondary lepton level, it would not 
be possible to predict the shape of these distributions accurately enough to determine M w . 
We refer the reader to Ref. [11] for a detailed discussion of the phenomenology of the W 
boson physics at the Tevatron collider. 

Similarly, at RHIC, a resummation calculation is needed to reliably predict the distribu- 
tion of the secondary leptons from the W boson decay. As we already mentioned, the reason 
is the necessity to deduce spin-dependent parton luminosities based on the observation of 
the secondary charged lepton only. For RHIC purposes, the resummation calculation at the 
lepton level has to be modified to include the spin dependence, which has been performed in 
Ref. [12]. In this paper, we will apply the calculation in Ref. [12] to model single-spin asym- 
metries of W boson production in realistic RHIC conditions. We put the emphasis on the 
discussion of the asymmetries Aiiye) and Aiiprt) for the distributions in the rapidity yt and 
transverse momentum prt of the observed charged lepton. We will argue that these directly 
observed asymmetries provide a viable alternative to the commonly discussed asymmetry 

Mvw)- 

The paper is organized as follows. In Section II, we briefly summarize the resummation 
formalism used for this study. A next-to-leading order (NLO) calculation is also presented 
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for comparison to the resummation calculation. Since good understanding of spin-averaged 
cross sections is needed to derive polarized PDFs from the spin asymmetries, Section III dis- 
cusses the uncertainties due to the imperfect knowledge of unpolarized parton distributions, 
as well as the potential for RHIC experiments to reduce these uncertainties. Section IV 
demonstrates that a correct model of multiple parton radiation is important to describe 
the differential distributions of the final-state leptons. We then employ the resummation 
calculation to realistically evaluate the feasibility of the reconstruction of the asymmetry 
Al(uw) from the observed charged lepton data. Section V shows that the spin asymmetries 
for the distributions of the observed secondary leptons provide an attractive alternative to 
the asymmetry Al(i/w)' they can be measured directly and discriminate efficiently between 
the different PDF sets. We also show that the RHIC experiments can explore the spin 
dependence of the nonperturbative contributions from the multiple parton radiation. The 
main findings of this paper are recapped in the conclusion. 

II. THEORETICAL ESSENTIALS 

The primary goal of this paper is to discuss production of W ± bosons, which will be used 
at RHIC to obtain information about the longitudinally polarized parton distribution func- 
tions (PDFs) in the quark sector. Due to the unique feature of the maximal parity violation 
in the Wqq coupling, the parton-level cross sections for spin-dependent W boson production 
have pronounced single-spin asymmetries. This feature gives W boson production an edge 
on parity-conserving processes, in which non-trivial hadronic dynamics can be probed only 
through the measurement of more complex double-spin asymmetries. 

Denote the unpolarized distribution of a parton a in a nucleon A as f a /A{£, A*f) and the 
longitudinally polarized distribution of a in the nucleon A as Af a /A(£, /Af)- Then, 



where fh a /h A (£,a,HF) is a helicity-dependent parton distribution function, i.e., a probability 
of finding a parton a with the momentum = £ a PA an ^ helicity h a in a hadron A with 
the momentum p A and helicity Ha- The PDFs depend on the factorization scale fj,p, which 
in our calculation is assumed to coincide with the QCD renormalization scale. To find the 



fa/A^,^F) = /+/+(£, Hf) + f -/+{£, Hf), 

A/„/a(C,M = /+/+(£> M - f -/+(£, hp), 



(2) 



(3) 
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polarized PDFs from an experiment with one longitudinally polarized beam, one measures 

the unpolarized cross section, 

da _ 1 (da(p^p -> WX) da(p^p -> WX)\ 

Ti~2 \ di + d~i ) ' (4) 



as well as the single-spin cross section, 

- n _, W X) \ 

(5) 



dA L a _ 1 (do(p^p -> WX) dcr(p^p -> WX) \ 



Here £ denotes any kinematic variable, such as the rapidity yw of the W boson. The single- 
spin asymmetry is defined as 

M ® = -we- (6) 

In Ref. [12], we presented fully differential unpolarized cross sections da/d^, single-spin 
cross sections dALd/d!;, and double-spin cross sections dA^ajd^ for production and de- 
cay of 7*,W ± , and Z° bosons. To reliably predict the rate at any point of the momentum 
phase space, these cross sections include the all-order sum of large soft and collinear loga- 
rithms, which dominate at small transverse momentum qx of the vector boson. This all-order 
sum is combined, without double counting, with the finite-order 0(ots) cross section, which 
dominates at large qr- The single-spin resummed cross sections in the narrow width ap- 
proximation were first obtained in Ref. [13]. The main result of Ref. [12] is the (more 
involved) derivation of the resummed cross sections that also accounts for the decay of the 
vector bosons and angular distributions of the final-state leptons. Our result was derived in 



the MS factorization scheme [14-16]; therefore, it is fully compatible with the MS parton 
distributions. 

The resummed cross sections from Ref. [12] are implemented in the numerical resumma- 
tion program Legacy and Monte-Carlo integration program RhicBos, 1 which were used 
to produce the results discussed in this paper. Legacy is a C++/Fortran program that 
quickly and accurately generates the resummed cross sections on a grid of points in the 
kinematical phase space. RhicBos reads in the cross section grids, performs their multi- 
dimensional integration, and produces the output that is optimized for RHIC specifics. Since 
the finite-order cross section is characterized by a smaller theoretical uncertainty at qr > Q, 



1 The Fortran code and input grids for RhicBos can be downloaded from 
http: / /hep. pa. msu.edu/~nadolsky/RhicBos/. 
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RhicBos switches from the resummed cross section (W^ + F) to the finite-order cross section 
at the point where the resummation cross section becomes smaller than the finite-order cross 
section at qr above the position of the maximum in da/dqr- To compare the resummation 
calculation to the finite-order approach, we have also calculated the finite-order O(as) cross 
section (next-to-leading order, or NLO, cross section) with the help of a phase space slicing 
method. This method introduces a separation transverse momentum scale q^ p , below which 
the integral of da jdq\ over the region < q\ < (<2^ ep ) 2 is calculated analytically, using the 
small-gr approximation in Eq. (78) of Ref. [12]. The integral over the region q\ > (q'^ p ) 2 
is integrated numerically using the exact 0(as) cross section. As we will show later, in the 
NLO approach, predictions for several important distributions depend strongly on the value 
of the auxiliary parameter q^ p , so that the resummed cross sections must be used in order 
to model those distributions. Some other distributions, such as da/dQ 2 , are not sensitive 
to q^ p (up to higher-order corrections), so that the resummation and NLO calculations give 
close results for such distributions. 

The electroweak parameters in this paper were obtained in the on-shell scheme according 
to the procedure described in Section II. B of Ref. [12] for the Fermi constant Gf = 1.16639 x 
1(T 5 GeV" 2 , W boson mass M w = 80.419 GeV, and Z° boson mass M z = 91.187 GeV. 
The width IV of the W boson was evaluated as T\y = 3G F M^/(2irV2). The finite-order 
and asymptotic cross sections were calculated for the QCD factorization scale \ip = Q. 
The resummed cross section (see Eqs. (79) and (80) in Ref. [12]) was calculated for the 
canonical values of the scale parameters, C\ = C3 = 2e~ 7E , where ~ 0.577 is the Euler 
constant, and Ci — \. Unless stated otherwise, the nonperturbative Sudakov factor in the 
resummed cross sections was evaluated by using the recent Gaussian parametrization [17] 
derived from the global analysis of transverse momentum distributions in unpolarized vector 
boson production. The projected statistical errors SA L in the measurement of the single-spin 
asymmetries were estimated according to Eq. (13) in Ref. [1]: 



where N = Ca is the number of the events for the given luminosity C and cross section a. 
Unless stated otherwise, the assumed integrated luminosity is C = 800 pb -1 , and polariza- 
tion P of the beam A is 70%. 




(7) 
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III. PHYSICS POTENTIAL IN UNPOLARIZED pp COLLISIONS 



To extract the single-spin cross section dAj^a/dC, from the single-spin asymmetry Al(£), 
the spin-averaged cross sections da/dt; must be known well. Obviously, the uncertainty in 
the knowledge of da/d^ must be small enough as compared to the targeted uncertainty 
in the measurement of dA^a/d^. Therefore, an interesting question is: how well do we 
currently know the unpolarized cross sections at RHIC energies, and how important is their 
measurement in the upcoming RHIC experiments? The goal of this section is to argue 
that at present the uncertainties in the unpolarized cross sections at RHIC are sizeable; 
to reliably predict the unpolarized rate, these uncertainties must be reduced by dedicated 
RHIC measurements. 

Table I lists the unpolarized cross sections and event rates for the processes pp — > (W ± — > 
tv()X and pp — > (Z° — > £l)X at RHIC for various values of center-of-mass energies y/s 
and integrated luminosities C. The cross sections were derived at NLO using the updated 
CTEQ6 parton distribution functions [18, 19]. The cross sections a are calculated using the 
best-fit PDF parametrization CTEQ6M. The uncertainty 5<jpdf in the cross section a is 
generated by the uncertainty in the unpolarized parton luminosities, which are themselves 
phenomenological functions known with a finite precision. This uncertainty is evaluated 
with the help of the Hessian matrix analysis employed in Ref. [18]. Specifically, 5<jpdf 
corresponds to the maximal variation of a for all possible PDF parametrizations lying in the 
acceptable range of \ 2 ■ According to Table I, the relative PDF uncertainties are of order 
25% at y/s = 200 GeV and 7 — 10% at y/s = 500 GeV. These are the errors that exist prior 
to the measurements at RHIC. 

The above PDF uncertainties are large because W and Z boson production at the pp 
collider RHIC probes the sea quark distributions (primarily u(x) and d(x)) at x > 0.1. For 
instance, the major part of W boson events occurs in the central-rapidity region yw ~ 
and corresponds to the Born-level momentum fractions x = Q/y/s close to 0.4 and 0.16 
at y/s = 200 GeV and y/s = 500 GeV, respectively. At such x the sea quark PDFs are 
not as constrained by the existing (mostly DIS) data, which better determine the valence- 
dominated distributions u(x) and d(x). Therefore, RHIC provides the information that is 
complementary to the DIS experiments at HERA and collider experiments at the pp collider 
Tevatron. 
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This point is illustrated by Fig. 1, which shows the relative PDF uncertainty SappF/v 
in W ± and Z° boson production at various existing colliders. To compare the potential 
of pp and pp colliders of the same center-of-mass energy, this figure also shows the PDF 
uncertainties at fictitious pp colliders with y/s = 200 and 500 GeV. According to the figure, 
the PDF uncertainties at RHIC substantially exceed the uncertainties at the Tevatron and 
LHC. They are also up to 2.5 times larger than the PDF uncertainties at the fictitious pp 
colliders of RHIC energies. The PDF uncertainty is the largest in W~ boson production at 
RHIC, which is sensitive primarily to the du combination of the parton distributions. 

On the other hand, for the planned integrated luminosities, RHIC will measure the spin- 
averaged rates with better accuracy, particularly in W + boson production. At t/s = 500 
GeV and C = 800 pb -1 , the 2a relative statistical errors {2N~ 1 / 2 ) in the rate are 0.6%, 1.2%, 
and 2% in W + , W~ , and Z° boson production, respectively. At y/s = 200 GeV and C = 
320 pb -1 , the respective errors are 10%, 18%, and 44%. We see that the projected statistical 
errors in W boson production are smaller than the current PDF errors. Therefore, we 
conclude that the measurement of the spin-averaged cross sections at RHIC will drastically 
reduce the uncertainty in the unpolarized parton distributions in the probed region of x. We 
emphasize that the measurement of the spin-averaged cross sections should be an important 
part of the RHIC physics program, since it complements information about the unpolarized 
PDFs from the other colliders and reduces systematic uncertainties in the measured polarized 
cross sections. 

IV. RECONSTRUCTION OF SPIN ASYMMETRIES FROM RHIC DATA 
A. Impact of multiple parton radiation on differential distributions 

To extract useful information from the raw data in W boson production, it is necessary to 
unfold the kinematical acceptance of the signal events, which can only be reliably predicted 
by a resummation calculation, similar to what has been done for the W boson physics at 
the Tevatron [11]. The correct lepton-level calculation is even more indispensable given that 
RHIC cannot register the transverse energy carried by the decay neutrino. 

For definiteness, let us concentrate on the production of W + bosons. Of course, all 
conclusions in this section hold for W~ boson production as well. Fig. 2 shows the transverse 
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momentum distributions of the W + boson and charged lepton from its decay, as predicted 
by the resummation package Legacy-RhicBos. Here and after, the numerical results 
correspond to the cm. energy y/s = 500 GeV and integrated luminosity £ = 800 pb -1 . 
It is evident that most of the W bosons are produced with small, but non-zero transverse 
momenta. Such non-zero qx is acquired through radiation of soft and collinear partons, 
which cannot be approximated by finite-order perturbative calculations. In order to obtain 
reliable predictions for differential cross sections, dominant logarithmic terms otg ln m (q^/Q 2 ) 
(where < m < 2n— 1) associated with such radiation should be summed through all orders 
of the perturbative series. 

Furthermore, consider the distribution da/dpre in Fig. 2(b). If the transverse momen- 
tum qT and width Tw of the W boson were negligible (as in the Born-level calculation), 
da/dpre would show a singular Jacobian peak exactly at pxe = M w /2. The smearing of 
the Jacobian peak in Fig. 2(b) is due to the multiple soft gluon radiation and non-zero 
width of the W boson. The example of da/dpre illustrates that multiple parton radiation 
and heavy boson decay may drastically alter the Born-level predictions. The resummation 
effects are important if the distribution is sensitive to directly (as da/dpre), or if the 
phase space is constrained in a region with slower convergence of perturbation theory (as 
da/dyi at large \ye\). In the remainder of this section, we utilize this accurate calculation to 
evaluate feasibility of the measurement of the boson-level asymmetry A L (y w ), which is com- 
monly discussed in literature as the most convenient observable to probe polarized sea quark 
distributions. Many aspects of our discussion are based on the analysis of two-variable dis- 
tributions d 2 o /{dpridye) in the transverse momentum and rapidity of the observed charged 
lepton. These distributions for the unpolarized and single-spin W + and W~ boson produc- 
tion are shown in Fig. 3. 



B. W boson rapidity at the Born level 

Consider the Born-level expression for the following asymmetry in W + boson production 
induced by the up quark and down antiquark via ud — > W + — > 

A / _ dA L a/(dQ 2 dy w d cos 6*) 

da / (dQ 2 dywd cos 9*) 
-Au(x A )d(x B )(l + cosr) 2 + Ad(x A )u(x B )(l - cosr) 2 

u(x A )d(x B )(l + cos6>*) 2 + d(x A )u(x B )(l - cos6>*) 2 ' ' 
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Here xa,b = (Q / y/s)e ±yw , and 9* is the polar angle of the charged lepton in the rest frame 
of the W boson, with the z axis pointing in the moving direction of the polarized proton 
beam. The single-spin asymmetry Al(vw) m the rapidity of the W boson can be obtained 
from Eq. (8) by integrating cos#* out, i.e., integrating over the moving direction of t + and 

A , x I-i dQ?d d yt L d cos e* dcosd * -Au{x A )d{x B ) + Ad{x A )u{x B ) 

The resulting expression for Ai(yw) is commonly discussed in literature as an example of 
a simple observable, which is straightforwardly sensitive to Au(x) and Ad(x). Similarly, 
the Born-level asymmetry A L (y w ) in W~ boson production is straightforwardly sensitive to 
Au(x) and Ad(x). It has been also argued in Ref. [22] that 0(as) corrections cancel to a 
good degree between the numerator and denominator of Ai(yw)- Unfortunately, as we will 
now demonstrate, the measurement of Al(uw) is obstructed by the effects of W boson decay 
and the limited detector acceptance, which seriously complicate both the reconstruction of 
yw and determination of A L {y w ) from the observed data. 

Let us first discuss the determination of the rapidity of the W boson. Given that neither 
of the two RHIC detectors measures the missing energy, the four-momentum of the W boson 
cannot be deduced from the momenta of its decay products. Therefore the rapidity of the 
W boson cannot be measured directly. 

Despite the impossibility to measure yw in general, it has been proposed [1, 5] to statis- 
tically choose the correct yw in a certain kinematical region based on several assumptions 
about the QCD dynamics of the process. Let us outline the main idea of this method. 

Denote the rapidity and transverse momentum of the lepton in the lab frame as ye and 
pre, respectively. Similarly, denote the corresponding variables in the rest frame of the W 
boson as y[ and p' Ti , respectively. In the naive Born approximation, the W bosons have zero 
width (Q = M w ) and zero transverse momentum (qx = 0). Therefore, the lepton variables 
in the different reference frames are related as 

vt = y'e + yw, (io) 

Vti = qr+Pre = Pre, ( n ) 
T 

q T = 



10 



where y' £ and p' Ti are the functions of the polar angle 9* of the lepton in the W boson rest 
frame: 



The comments below the equations show the assumptions involved in obtaining the rightmost 
expressions in Eqs. (11) and (13). 

For a given set of yt and pre, Eqs. (10)-(13) imply two solutions y + and y_ (corresponding 
to 9* and n— 9*) for yw- One of them can be correctly chosen if the magnitude of the rapidity 
ye of the charged lepton is large (i.e., if the lepton is observed in the forward or backward 
rapidity region). In these special regions, one of the two solutions for yw can always be 
discarded, because it takes a value outside of the allowed rapidity range of the collider. In 
other words, the chosen solution is the one with the smaller magnitude. However, in reality, 
qx never vanishes, as predicted by the resummation calculation (see Fig. 2(a)). Similarly, 
Yw is not zero. Therefore, the above approximation never holds exactly, and the size of the 
error is determined by the difference between the Born-level and exact dynamics. 

By comparing the leading order, NLO, and resummation calculations, we have found 
that the approximation holds reasonably well when qx <C M w , p' Tl m p Til and yg is large. 
This point is illustrated in Fig. 4, which shows the statistically averaged differences ((\y± — 
yw\)&) I (&) between the exact yw and each of the two solutions y± in the two-dimensional 
plane of yi and pre, predicted by the resummation calculation. It can be seen that the 
difference is small for one of the solutions in the respective region yt > 1 or yi < — 1, and 
Pre < 25 — 30 GeV. At the same time, the difference is unacceptably large for the other 
solution. Correspondingly, in each special region the distribution in the correct solution 
y + or y_ closely reproduces the distribution in true yw- For example, Fig. 5 shows that 
the distributions da/dy_ agrees well with the distribution da/dyw in the region 10 GeV < 
pre < 30 GeV and 1.2<y t < 2.4. 

Except for the special regions in Fig. 4, the approximate solutions y± do not agree well 
with the exact yw- An important feature to note is that in the special regions the event rate 
is small. As shown in Fig. 3, most of the event rate comes from pre ~ and yn ~ 0. 

Furthermore, the magnitude of the error in yw depends on the shape of the cross section, 




1, 1+COS0* 

- In , 

2 l-cos#* 

Q ■ n* 

— sm9 = 
2 

Q = M w 



M w 



2 



sin 9*. 



(12) 



(13) 
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which, of course, is unknown a priori in the polarized case. 
C. Distortion effects on Al(2/w) 

Let us now turn to another problem in the measurement of Al(i/w)' distortions in the 
deduced shape of Al(vw) due the experimental cuts. To obtain Al(i/w), the experimental 
cross sections have to be integrated over all possible directions of the lepton motion, i.e., over 
cos 9* in the interval —1 < cos 9* < 1 (see Eq. (9) for an example of such integration in the 
Born-level analysis). Unfortunately, the acceptance of RHIC detectors does not permit such 
integration, so that the resulting asymmetry has the residual dependence on the acceptance 
range in cos#* and does not agree with the true asymmetry Ai(yw). 

Indeed, cos9* is related to the rapidity of the lepton ye via Eqs. (10) and (12): 

1 i + cos e* 

2 1 — costr 

Therefore, the restrictions on the range of ye reduce the range of integration over cos#*. 
Specifically, the Phenix muon detectors are able to detect a muon if the muon rapidity 
y^ and azimuthal angle M satisfy 1.2 < \y^\ < 2.4 and < M < 2tc. The Phenix 
electromagnetic calorimeter can register electrons with \y e \ < 0.35 and < e < n. The 
Star detector registers electrons in the barrel electromagnetic calorimeter that covers \y e \ < 
1.0 and < <j) e < 2n. In addition, a selection cut pre > 10 — 20 GeV is expected to be 
imposed on the charged leptons to suppress background contributions, which also constrains 
sin#* through Eqs. (11) and (13). 

The effect of the above cuts is demonstrated in Fig. 6, which shows the single-spin asym- 
metry in W + boson production calculated without constraints on y^ and pti (solid line), as 
well as with the constraints 1.2 < \ye\ < 2.4, pre > 20 GeV (circles) and \yg\ < 1, pre > 20 
GeV (boxes). The results shown here are derived using the resummation calculation. 2 Ac- 
cording to the figure, there is a substantial difference between the asymmetries calculated 
with and without experimental cuts. This difference arises due to the different dependence of 

2 For the asymmetries Al(dw), the resummation and NLO calculations give close predictions (within one 
estimated statistical error SAl), even though these calculations predict quite different shapes for the 
unpolarized and polarized differential cross sections. 
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the unpolarized and polarized cross sections on angular distributions of the charged leptons, 
which affects the asymmetry because cos^* is not integrated out completely. 

The findings in this section can be summarized as follows. The extraction of Aiiyw) from 
RHIC data is not an easy task, because (a) the direct measurement of yw is impossible, 
and (b) the reconstructed Ai(yw) is distorted by the experimental acceptance cuts. In 
particular, the two-fold ambiguity in the determination of yw is unavoidable at ye ~ and 
PTe ~ Mw/2, i.e., in the region with the largest cross section. To take the full advantage 
of the large event rate at RHIC, we ought to seek an alternative to A^yw), which is free 
of the complications discussed above. In the next section, we turn our attention to the spin 
asymmetries of the lepton-level cross sections da /dye and da/dpxe, which may serve as such 
alternative observables. Based on the resummation calculation, we analyze the sensitivity 
of the lepton-level asymmetries to the choice of the experimental cuts. We also present 
theoretical predictions for various PDF sets. 

V. SINGLE-SPIN ASYMMETRIES AT THE LEPTON LEVEL 

A. Impact of the cuts and soft radiation on lepton-level cross sections 

Let us first separately discuss the unpolarized and polarized cross sections. To study the 
detector effects on the shape of various distributions, we will compare the cross sections for 
the kinematical acceptance of Phenix and Star detectors, as well as cross sections without 
any kinematical cuts. We also compare the resummation calculation to the NLO calculation, 
with the latter performed in the phase space slicing method for the separation parameter 
q^ p = 1.6 GeV. We concentrate on the distributions da/dpre, where the differences between 
the resummed and NLO predictions are the most visible. 

Fig. 7 presents the distribution da /pre for the positively charged leptons from the decay 
of W + bosons. The shape of the NLO cross sections (dashed curves) for pre around and 
above M w /2 (i.e., the Jacobian peak) is strongly affected by the arbitrary parameter q^ v 
employed in this calculation. Hence, the only reliable predictions for this part of the phase 
space are given by the resummed cross sections (shown as the solid curves). For smaller pre, 
the resummation calculation agrees with the NLO calculation, because it is formulated to 
match the NLO calculation in the region where the multiple gluon emission is not important. 
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Furthermore, without any kinematic cuts, the integrated rate of the resummation calculation 
agrees with that of the NLO calculation to a percent level. As shown in Fig. 7(b), the 
selection of the leptons in the large rapidity region distorts the shape of the pre distribution 
by suppressing the contributions near the Jacobian peak. On the other hand, the selection 
of the central-rapidity region \y e \ < 1.0 (cf. Fig. 7(c)) suppresses the contributions with 
Pre ;$ 10 GeV. Both features can be understood from the shape of the differential distribution 
d 2 a/ (dpTidye) given in Fig. 3(a). According to this figure, the small (large) rapidity \yg\ 
statistically corresponds to the large (small) transverse momentum pre- 

As a pp collider, RHIC produces different numbers of W~ and W + bosons because of 
different parton luminosities. Due to the different rapidity distributions for W + and W~ 
bosons, the rapidity cuts have a different effect on da/dpre in these processes. Fig. 8 shows 
the distributions da/dpTe in the process pp — > {W~ — > t~V()X. We find that the shape of the 
distributions without cuts shown in Fig. 8(a) is very close to the shape of the corresponding 
distribution in Fig. 7(a), even though the overall normalizations are obviously different. On 
the other hand, the shapes do not coincide in the presence of the cuts, in particular in the 
case of the Phenix selection cuts shown in Figs. 7(b) and 8(b). 

The pti distributions in the single-spin W + boson production are shown in Fig. 9. The 
cross sections dA L a/dpxe in W + boson production are negative, because the charged weak 
coupling is purely left-handed, and the dominant polarized PDF Au(x) is positive in the 
RHIC range of x. Again, we see that the NLO calculation cannot reliably predict the 
distribution of pre around and above Mw/2. For the Phenix kinematics, cf. Fig. 9(b), 
a sharper peak is developed in the distribution of pre around as compared to the 

unpolarized case, cf. Fig. 7(b). The different shape is caused by the different rapidity 
dependence in the unpolarized and single-spin cases, as shown in Figs. 3(a) and 3(c). Finally, 
the distribution da/dp T g for the charged leptons t~ in the W~ events is shown in Fig. 10. 

For the rapidity distributions da /dye without imposed cuts on pti, the NLO and re- 
summed cross sections are close to one another, in accordance with the cancellation of 
soft contributions in sufficiently inclusive observables. However, according to the discussion 
above, the resummation and NLO predictions for da /dye can be very different if cuts on pxe 
are imposed. 
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B. Single-spin asymmetries 



We now turn our attention to the single-spin asymmetry. Since the only straightforward 
signature of the W bosons at RHIC is the observation of the secondary charged leptons, 
it is important to understand these asymmetries at the lepton level. In Figs. 11 and 12, 
we show the lepton-level asymmetries Ai(ye) in the W + and W~ events, respectively, for 
various cuts on pre- In these figures, the positive rapidity direction is in the moving direction 
of the polarized proton beam. Similarly, Figs. 13 and 14 show the asymmetries A L {p T e) for 
various cuts on ye. All asymmetries are derived using the Gehrmann-Stirling PDF sets A 
and B [23] and GRSV-2000 valence-like PDF set [21] in the resummation calculation. We 
found by experimenting with other available PDF sets in Refs. [21, 23, 24] that the shown 
graphs characterize well the possible span in the asymmetry due to the different polarized 
PDF parametrizations. The projected statistical errors are estimated according to Eq. (7) 
for the bin sizes shown in the figure. 

According to the figures, the lepton-level asymmetries can be accurately measured for 
both W + and W~ bosons. These directly observed asymmetries can efficiently discriminate 
between different PDF sets; hence, they provide a viable alternative to the less accessible 
asymmetry A L (y w ). For both jy + and W~ boson samples (cf. Figs. 11 and 12), the power 
to discriminate between different PDF sets is larger in the region of the positive rapidity 
ye, especially for the W~ boson data. It is interesting to observe that, in different ranges 
of ye, the largest variations due to the choice of the PDFs appear in different ranges of pre- 
In W + boson production, the largest variation in the asymmetries A^{j>ti) is concentrated 
a t Pre > 30 — 40 GeV in the backward rapidity region ye < (Fig. 13(a) and (b)) and at 
Pre ^ 30 — 40 GeV in the forward rapidity region ye > (Fig. 13(c) and 13(d)). In W~ 
boson production the largest variation is observed for very forward leptons, 1.2 < ye < 2.4 
(Fig. 14(d)). If an additional condition pre > 20 GeV is imposed, other regions of ye 
(cf. Fig. 14(a)-14(c)) show a smaller sensitivity to different PDF sets. 

During the first years of RHIC spin program, it may not be possible to distinguish between 
the positively and negatively charged leptons from the W boson decay, particularly before 
the solenoid field of the Star detector is fully operational. It is therefore interesting to 
examine the possibility of measuring the single-spin asymmetries in the combined sample of 
W + and W~ bosons, without distinguishing between charged leptons and anti-leptons. Such 
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asymmetries are shown in Figs. 15 and 16. Since there are about three W + bosons produced 
per each W~ boson, these asymmetries are close to the W + boson asymmetries shown in 
Figs. 11 and 12. Therefore, these asymmetries will be able to distinguish between different 
PDF sets for a sufficiently high number of events. 3 However, they cannot completely replace 
the single-spin asymmetries for independent W + and W~ boson samples, which have better 
discriminating power and allow us to distinguish between contributions of different parton 
flavors. 

C. Nonperturbative contribution 

As was discussed above, the transverse momentum distributions for vector bosons are af- 
fected by the soft and collinear radiation, which can be described only by means of all-order 
resummation. In addition, the distributions at very small qx are sensitive to nonpertur- 
bative contributions characterized by large impact parameters b > 1 GeV -1 . As a result, 
the shape of the lepton-level distribution da/dpre. around its peak at about pre = Mw/2 
is affected by both perturbative and nonperturbative QCD radiation, in addition to the 
nonzero width of the W boson. While the perturbative soft and collinear contributions can 
be calculated order-by-order in the resummation formalism, the nonperturbative contribu- 
tions (associated with the nonperturbative transverse motion of partons inside the proton) 
can only be approximated by phenomenological parametrizations. In particular, based on 
the factorization properties of the resummed cross sections, one can hypothesize a simple 
relationship between the nonperturbative Sudakov factors in the unpolarized, single-spin, 
and double-spin vector boson production: 

^ ^ p\ single— spin 2 ( ^NP \ unpolarized ^ ^ p\ double— spin^j ' (l^) 

Furthermore, Snp is usually assumed to be independent of the type of the electroweak hard 
probe, so that it is the same in 7*, W, and Z boson production. If the relationship (14) for the 
spin dependence of Snp indeed holds in nature, the nonperturbative contributions in single- 
spin W boson production can be determined unambiguously once the parametrizations of 

3 Of course, the projected error bars in Figs. 15 and 16 (evaluated for C — 800 pb" 1 ) will have to be 
adjusted to agree with the actual integrated luminosity at the time of the measurement. 
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Snp on the right-hand side of Eq. (14) are available, e.g., from the global analysis of the 
distributions da/dqr in the unpolarized and double-spin 7* or Z° production. 

It is well known from the unpolarized studies (see, e.g., Refs. [10, 25, 26]) that the 
importance of the nonperturbative contributions in the massive boson production at high- 
energy colliders is strongly reduced as compared to Drell-Yan pair production in fixed-target 
experiments. Nonetheless, this dependence is not entirely negligible at small transverse mo- 
menta of W bosons, roughly for qx ~ 5 GeV or less. To give an idea about the typical 
size of the nonperturbative contributions, Fig. 17 shows the event rate for the single-spin 
cross section dAj^a /dpre at C = 800 pb" 1 , integrated over the Star acceptance range 
— 1 < He < 1- This rate was calculated using two parameterizations [17, 27] of the nonper- 
turbative part, which were found from the analysis of unpolarized vector boson production. 
Correspondingly, the resummed cross sections were calculated with the assumption that 
SNp\ single _ spin = S NP \ unpolarized . It can be seen that the sensitivity to the nonperturba- 
tive input is small, but, nonetheless, visible near the Jacobian peak. For comparison, we 
also included the NLO cross section calculated using the phase space slicing method. The 
NLO curve substantially deviates from the resummation curves, and, moreover, its shape 
can be drastically modified by varying the phase space slicing parameter q^ p ■ In contrast, 
practically all features of the resummation curve are determined unambiguously by the 
perturbative contribution, and the remaining small dependence on the nonperturbative con- 
tributions is pinned down by utilizing a phenomenological parameterization discussed above. 
Needless to say, RHIC should explore the spin properties of nonperturbative contributions 
by studying distributions da/dqr in the double-spin 7* and Z° production and da/dpre. in 
single-spin and double-spin W boson production. 



VI. CONCLUSION 



Due to the maximal parity violation in the weak current coupling, the longitudinal single- 
spin asymmetry of W boson production is strongly sensitive to the polarized parton distri- 
butions. At the center-of-mass energy of 500 GeV, W boson production events at RHIC can 
clearly discriminate between various polarized quark luminosities. As a pp collider, RHIC 
probes the sea quark PDFs in the region of relatively high x, where such PDFs are not ade- 
quately constrained even in the unpolarized case. Therefore, RHIC experiments also have a 
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good potential to reduce the uncertainties in the unpolarized sea quark PDFs in a dedicated 
study of spin-averaged cross sections. 

Since Phenix and Star detectors do not determine the missing transverse energy as- 
sociated with the neutrino from W boson decay, and because the charged leptons from the 
decay can only be observed in a part of the 4n solid angle, the reconstruction of the asym- 
metries at the level of W bosons is not easy to realize. As an alternative to the asymmetry 
Al(uw) for the W boson rapidity distribution da/dyw, we advocate the measurement of the 
asymmetries Al(z^) and Al(pt£) for the distributions da /dye and da/dpre in the rapidity 
ye and transverse momentum pre of the secondary charged lepton. These lepton-level asym- 
metries can be directly observed at RHIC detectors, they are sizeable, and they can be used 
to determine the polarized PDFs with high precision. 

To reliably predict the distributions of the charged leptons in all available phase space, a 
theoretical calculation is needed for all-order summation of large logarithmic contributions 
corresponding to W bosons with a small transverse momentum. At the lepton level, the 
largest effect of the logarithm resummation is observed in the distributions da/dpre and 
df±.L<J I dpTt in the proximity of the Jacobian peak (pre ~ M^/2). The importance of the 
resummation can be verified by observing a large difference between the resummed and 
NLO cross sections near the Jacobian peak, which happens only when the effect of soft 
gluon radiation is important. 

In this paper, we apply the state-of-art resummation calculation [12] for spin-dependent 
vector boson production to predict reliably the distributions of the observed leptons from W 
boson decay. In addition to the spin-averaged and single-spin cross sections, we presented a 
detailed study of single-spin asymmetries for the decay lepton in the presence of kinemat- 
ical cuts imposed by Phenix and Star detectors. Furthermore, we propose to study the 
distribution dA L a/dpTe at pxe ~ M w /2 to learn about the spin (in) dependence of nonper- 
turbative contributions in the CSS resummation formalism. The study of such contributions 
for different beam polarizations and bosons of different types ( r y*,W ± , Z°) will shed light on 
basic properties of QCD factorization and provide important clues about the nature of the 
intrinsic transverse motion of partons inside the proton. 



18 



Acknowledgments 

Authors would like to thank C. Balazs, D. Boer, G. Bunce, M. Grosse Perdekamp, 
S. Gupta, J. Kiryluk, N. Saito, M. Stratmann, W. Vogelsang, and members of the CTEQ 
collaboration for helpful discussions. We are grateful to F. Olness for his comments on the 
manuscript. We thank the organizers of RHIC Spin workshops, where preliminary results of 
this work were presented. The work of P. M. N. has been supported by the U.S. Department 
of Energy and Lightner-Sams Foundation. The research of C. P. Y. has been supported by 
the National Science Foundation under grant PHY-0100677. 



[1] G. Bunce, N. Saito, J. Soffer, and W. Vogelsang, Ann. Rev. Nucl. Part. Sci. 50, 525 (2000). 

[2] C. Bourrely and J. Soffer, Nucl. Phys. B423, 329 (1994). 

[3] C. Bourrely and J. Soffer, Nucl. Phys. B445, 341 (1995). 

[4] P. M. Nadolsky (1995), hep-ph/9503419. 

[5] L. C. Bland (STAR Collaboration) (2000), hep-ex/0002061. 

[6] D. P. Morrison et al. (PHENIX Collaboration), Nucl. Phys. A638, 565 (1998). 

[7] J. W. Harris (STAR Collaboration), Nucl. Phys. A566, 277c (1994). 

[8] T. Affolder et al. (CDF Collaboration), Phys. Rev. D64, 052001 (2001). 

[9] V. M. Abazov et al. (DO Collaboration), Phys. Rev. D66, 012001 (2002). 

[10] J. C. Collins, D. E. Soper, and G. Sterman, Nucl. Phys. B250, 199 (1985). 

[11] C. Balazs and C.-P. Yuan, Phys. Rev. D56, 5558 (1997). 

[12] P. M. Nadolsky and C.-P. Yuan (2003), hep-ph/0304001. 

[13] A. Weber, Nucl. Phys. B403, 545 (1993). 

[14] R. Mertig and W. L. van Neerven, Z. Phys. C70, 637 (1996). 

[15] W. Vogelsang, Phys. Rev. D54, 2023 (1996). 

[16] W. Vogelsang, Nucl. Phys. B475, 47 (1996). 

[17] F. Landry, R. Brock, P. M. Nadolsky, and C. P. Yuan (2002), hep-ph/0212159. 

[18] J. Pumplin et al., JHEP 07, 012 (2002), hep-ph/0201195. 

[19] D. Stump et al. (2003), hep-ph/0303013. 

[20] H. L. Lai et al. (CTEQ Collaboration), Eur. Phys. J. C12, 375 (2000). 



19 



[21] M. Gluck, E. Reya, M. Stratmann, and W. Vogelsang, Phys. Rev. D63, 094005 (2001). 

[22] T. Gehrmann, Nucl. Phys. B534, 21 (1998). 

[23] T. Gehrmann and W. J. Stirling, Phys. Rev. D53, 6100 (1996). 

[24] M. Gluck, E. Reya, M. Stratmann, and W. Vogelsang, Phys. Rev. D53, 4775 (1996). 

[25] C. T. H. Davies, B. R. Webber, and W. J. Stirling, Nucl. Phys. B256, 413 (1985). 

[26] P. B. Arnold and R. P. Kauffman, Nucl. Phys. B349, 381 (1991). 

[27] G. A. Ladinsky and C.-P. Yuan, Phys. Rev. D50, 4239 (1994). 



20 



Table I: The unpolarized cross sections (in pb) and numbers of events N for the processes pp — ► 
(W^ — ► lv()X and pp — > (Z° — > at RHIC. Here i and f£ are leptons of one lepton generation, 
and the decay branching ratios were calculated to be Br (W ± -> tv t ) = 0.11, Br (Z° -> it) = 
0.034. The numbers in parentheses show the relative uncertainties. Also shown are the Born-level 
momentum fractions x = Qj^fs corresponding to yw = (or yz = 0). The uncertainties dapoF 
in the cross sections a are due to the uncertainties in the parton distribution functions estimated 
according to the CTEQ6 analysis [18, 19]. The uncertainties in N are purely statistical. 
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Figure 1: Uncertainties due to the PDF errors in the unpolarized cross sections of W and Z boson 
production at various colliders. The shown uncertainties are for RHIC at yfs = 200, 500, and 650 
GeV; Tevatron Run-2; LHC; and fictitious pp colliders at = 200 and 500 GeV. 
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q T ,GeV p Tl ,GeV 
(a) (b) 

Figure 2: The unpolarized cross section for the process pp — > (W + — > ^ + z^)X at = 500 GeV as 
a function of (a) the transverse momentum qr of the W + boson and (b) the transverse momentum 
Pti of the charged lepton, as predicted by the resummation calculation. CTEQ5M PDFs [20] were 
used. 
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(c) A L pp -» (W+ -» (d) A iP p - (w- - r 

Figure 3: Resummed two- variable distributions d 2 a/dpTedye for (a), (b) unpolarized and (c), (d) 

single-spin cross sections in W + and W~ boson production. The unpolarized and single-spin cross 

sections are calculated using CTEQ5M [20] and GRSV-2000 [21] standard PDF sets, respectively. 

A high-resolution version of this figure is included in the online copy of the preprint available at 

http:/ /hep.pa.msu.edu/~nadolsky/RhicBos/home.html#RhicBos_references. 
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Figure 4: The average differences between the exact value of yw and two approximate solutions 
y± of Eqs. (10)-(13) as a function of the charged lepton's rapidity yi and transverse momentum 
pxe- The resummation calculation and CTEQ5M PDFs [20] were used. The numbers inside the 
"butterfly" shapes show the magnitude of the differences \y + — yw\ and |y_ — yw\- According 
to the figure, the solutions y + and y_ provide good approximations in the limits y£ <C and 
yi S> 0, respectively The dashed lines approximately mark the regions of the shown rapidity 
differences. A more detailed color version of the figure is available in the online copy of the preprint 
at http: / /hep.pa.msu.edu/~nadolsky/RhicBos/home.html^RhicBos_references. 
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Figure 5: Comparison of the distributions da/dy- and da/dyw in the region of the validity of the 
approximate solution y_, as predicted by the resummation calculation. The constraints 1.2 < yi < 
2.4 and pri > 10 GeV imitate kinematical cuts for a measurement with the Phenix muon detector. 
The cut pti < 30 GeV suppresses contributions with large deviations of y_ from yw (cf. Fig. 4(b)). 
CTEQ5M PDFs [20] were used. 
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Figure 6: Dependence of the asymmetry Al(i/w) on the cuts imposed on the momentum of the 
observed charged lepton in the process A^pp — > (W + — ► l + ve)X, as predicted by the resummation 
calculation. The GRSV-2000 standard set [21] of the polarized PDFs was used. The projected 
statistical errors in Al{vw) (shown by error bars) are estimated with the help of Eq. (7) for the 
bin sizes shown in the figure. 
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Figure 7: Unpolarized resummed (solid) and NLO (dashed) cross sections da/dp^g in the process 
pp — ► (W + — > ^ + ^)X for various cuts on the lepton's rapidity. CTEQ5M PDFs [20] were used. 
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Figure 8: Unpolarized resummed (solid) and NLO (dashed) cross sections da/dp^g in the process 
pp — ► (W~ — > t~i>f)X for various cuts on the lepton's rapidity. CTEQ5M PDFs [20] were used. 
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Figure 9: Single-spin resummed (solid) and NLO (dashed) cross sections dA^a/dpTe in the process 
Alpp — > {W + — > £ + vi)X for various cuts on the lepton's rapidity. Gehrmann-Stirling PDF set A 
[23] were used. 
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Figure 10: Single-spin resummed (solid) and NLO (dashed) cross sections dA^a / dpxe in the process 
Alpp — > — > £~i>i)X for various cuts on the lepton's rapidity. Gehrmann-Stirling PDF set A 
[23] was used. 
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Figure 11: Asymmetries Al(v{) for various selection cuts on pxe in W + boson production, as 
predicted by the resummation calculation. The asymmetry is shown for the Gehrmann-Stirling 
PDF sets A (solid) and B (dashed), as well as for the GRSV valence-like PDF set (dotted). 
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Figure 12: Asymmetries Ai,{yg) for various selection cuts on pti in W~ boson production, as 
predicted by the resummation calculation. The asymmetry is shown for the Gehrmann-Stirling 
PDF sets A (solid) and B (dashed), as well as for the GRSV valence-like PDF set (dotted). 
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Figure 13: Asymmetries Ai,(pTe) for various selection cuts on yi in W + boson production, as 
predicted by the resummation calculation. The asymmetry is shown for the Gehrmann-Stirling 
PDF sets A (solid) and B (dashed), as well as for the GRSV valence-like PDF set (dotted). 
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Figure 14: Asymmetries Al(pti) for various selection cuts on yg in W~ boson production, as 
predicted by the resummation calculation. The asymmetry is shown for the Gehrmann-Stirling 
PDF sets A (solid) and B (dashed), as well as for the GRSV valence-like PDF set (dotted). 
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Figure 15: Asymmetries A^d/i) for various selection cuts on pri for the combined sample of W + 
and W~ bosons, as predicted by the resummation calculation. The asymmetry is shown for the 
Gehrmann- Stirling PDF sets A (solid) and B (dashed), as well as for the GRSV valence-like PDF 
set (dotted). 
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Figure 16: Asymmetries A^ipxe) for various selection cuts on yi for the combined sample of W + 
and W~ bosons, as predicted by the resummation calculation. The asymmetry is shown for the 
Gehrmann- Stirling PDF sets A (solid) and B (dashed), as well as for the GRSV valence-like PDF 
set (dotted). 
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Figure 17: The single-spin charged lepton transverse momentum distribution dALa /dpre for W + 
boson production discussed in the main text. The nonperturbative parts of the resummed cross 
sections were approximated by the Ladinsky-Yuan [27] (solid) and Brock-Landry-Nadolsky-Yuan 
[17] (dashed) phenomenological parameterizations from unpolarized vector boson production. The 
NLO cross section is shown for comparison as a dotted line. The distributions are calculated using 
the Gehrmann- Stirling PDF set A [23]. 
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